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Abstract 

We examine the modulation effects of thicknesses of both a photoactive layer (a bulk-heterojunction (BHJ) of 
poly(3-hexylthiophene) and [6,6]-phenyl-C61 -butyric acid methyl ester (P3HT:PCBM)) and an optical spacer of a 
transparent metal oxide, for power conversion efficiency optimization of organic photovoltaic devices. The 
redistribution of the optical intensity at the photoactive layer via the thickness modulation of both layers is taken 
into account, to produce three-dimensional (3D) plots as a function of both layer thicknesses of 0 to 400 nm range 
(5 nm step), for the device efficiency optimization. The modulation pattern of absorption is produced in the 3D plot 
as scanning the thicknesses of both layers as a result of modulation of interference between incoming and reflected 
light, which can be secured by changing the effective optical path length between two electrodes of a photovoltaic 
device. It is also seen that the case of inserting the spacer of the higher refractive index demands finer adjustment 
of the spacer layer thickness to achieve the optimum device efficiency. 

In addition, the series resistance of the photoactive layer of the thickness range of 0 to 70 nm is taken into account 
to provide the 3D plots as a function of the scanned thicknesses of both layers. Inclusion of the series resistance of 
the photoactive layer, which is also the function of its thickness, in the simulation, indicates that the series 
resistance can influence qualitatively the dependence of power conversion efficiency (PCE) on the thicknesses of 
both layers. We also find that minimization of series resistance, e.g., by device annealing, allows not only the 
relevant voltage to increase but also the optimum thickness of the photoactive layer to increase, leading to more 
absorption of light. 
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Background 

The past several decades have witnessed the rapid devel- 
opment of an organic photovoltaic device. Enormous ef- 
forts have been devoted to the improvement of organic 
photovoltaic efficiencies [1-6] with the recent report of 
power conversion efficiency (PCE) of about 12% [7-11]. 
Particularly, it is known that bulk-heterojunction (BHD- 
based organic photovoltaic devices exhibit advantages 
such as the low cost, its flexibility, its light weight, and 
the solution processibility for its fabrication. 
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However, the limited electron-hole transport character- 
istics of a BHJ photoactive layer do not simply encourage 
the active layer elongation that can benefit from the 
expansion of absorption region, because the elongated 
transport distance increases the rate of recombination of 
photo-generated charge carriers [12]. To enhance the de- 
vice performance, various device structures that employ 
optical spacers [13,14], photonic crystals [15], and local- 
ized surface plasmon resonance are proposed. More 
interestingly, optical spacers inserted between two elec- 
trodes of a photovoltaic device have been studied to re- 
distribute light intensity profiles across the photoactive 
layers where absorption occurs to maximize the absorp- 
tion. This redistribution is attributed to optical interfer- 
ence changes, resulting from the changes in an optical 
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path length between two metal electrodes that collect 
photo-generated carriers. 

Previous studies on optical spacers, which show en- 
hancement of absorption at the photoactive layer via the 
insertion of optical spacers, however, have used a number 
of discrete values of spacer thicknesses over the limited 
range for given thicknesses of the photoactive layers. In 
addition, one may not simply estimate the device per- 
formance, i.e., short-circuit current as thicknesses of both 
photoactive and spacer layers are modulated due to the 
following two reasons: first, both thickness modulations 
(two different materials) complicate the redistribution of 
light intensity at the photoactive layer via the modulation 
of the effective optical path length; and second, the 
photoactive layer thickness also governs the series resist- 
ance of the device, which consists of both the active layer 
bulk resistance and its contact resistance with the elec- 
trode [16-19]. For example, the increase in the photo- 
active layer thickness that is expected to expand the 
absorption region may not simply result in the enhance- 
ment of the PCE due to the resistance change. 

In this paper, we theoretically investigate the performance 
efficiency of a BHJ organic photovoltaic film comprising a 
blend of semiconducting polymer (poly(3-hexylthiophene): 
P3HT) and fullerene derivative ([6,6]-phenyl-C61-butyric 
acid methyl ester: PCBM). This photoactive film is sand- 
wiched between the charge separating electrodes, i.e., the 
hole collecting electrode of poly(3,4-ethyl-enedioxythio- 
phene):poly(styrene-sulfonate): PEDOT:PSS-indium tin 
oxide (ITO) glass, and the electron collecting electrode 
of aluminum (Al). Insertion of a metal oxide optical spa- 
cer between the Al electrode and the photoactive layer 
(P3HT:PCBM) is considered to optimize the optical 
intensity distribution over the photoactive region. 

By scanning the thicknesses of both the photoactive 
and optical spacer layers from 0 to 400 nm with a 5 nm 
step, we provide a three-dimensional (3D) plot of the 
short-circuit current (J sc ) at the photoactive layer, taking 
into account the expansion of absorption region and light 
intensity redistribution at the layer. The modulation 
behavior of J sc is observed along the scanned thickness of 
each layer over the entire range scanned in the 3D graph 
as a result of optical field redistribution caused by the ef- 
fective path-length change between two electrodes. The 
modulation behavior becomes more faded out at larger 
thickness of the photoactive layer. Optical spacers of dif- 
ferent refractive index materials, i.e., ZnO and Ti0 2 , are 
also examined, revealing that the higher index spacer de- 
mands the finer adjustment of the spacer layer thickness 
for device efficiency optimization. 

The series resistance (R) which is known to vary with 
the photoactive layer thickness [16] is used to normalize 
Jsc to represent the device performance (y\ = J SC IR). Inclu- 
sion of such thickness-dependent resistance in the 3D 



plots of t] produces the qualitative change in the 3D plot 
shape of the normalized J sc , implying that the thickness- 
dependent resistance causes the qualitative change in 
the dependence of device performance on both thick- 
nesses of the active layer and the spacer. It is therefore 
seen that the reduction of series resistance permits not 
only the relevant voltage to increase but also the optimum 
thickness of the photoactive layer to increase leading to 
more absorption that takes place. 

Methods 

There are two commonly used methods to calculate distri- 
butions of optical fields in a multilayer structure of thin 
films where interference between incident and reflected 
light occurs. One of them is by using a transfer matrix 
method based on matrices of optical propagation during 
the constituent media and those of interface with field 
boundary conditions [20]. The other one is the finite dif- 
ference time domain (FDTD) method, which has recently 
been widely used for solving Maxwell's equations for a 
given optical system of a complex structure. 

A FDTD method relies on numerical calculation with 
discrete meshes in both space and time. A discrete mesh 
made of so-called Yee cells is used to describe the elec- 
tromagnetic fields and structural materials of interest. 
Maxwell's equations are solved discretely in time, where 
the time step used is dominated by a mesh size through 
the speed of light. We use the software called FDTD Solu- 
tions (Lumerical Solutions, Inc., Vancouver, BC, Canada) 
for calculation of optical field distribution within the 
photovoltaic cell of interest. This calculation includes the 
interference effects between optical fields propagating in 
the forward direction and those in the backward direction 
within the photovoltaic device. 

Scanning of thicknesses of device constituent layers 
changes a condition for optical interference between in- 
cident and reflected light, thus allowing the optimum 
condition to be achieved for optimum field distribution 
for energy harvest efficiency. Figure 1 shows schematic 
of the organic photovoltaic structure considered in 
the simulation, where we scan the thicknesses of both 
the active layer (P3HT:PCBM) and the optical spacer, 
while keeping constant the thicknesses of the other 
layers, i.e., an ITO and PEDOT:PSS layers as 180 and 
50 nm, respectively. 

At first, the number of photons absorbed in the 
photoactive layer of the BHJ photovoltaic device is com- 
puted with an aid of the software FDTD Solutions, while 
the set of required optical parameters of the materials 
such as the relevant absorption coefficient are obtained 
from [21], and those for the sunlight spectrum from 
ASTM [22]. As an optical spacer layer, a zinc oxide 
(ZnO) which is a transparent electron acceptor with the 
refractive index of about 1.4 at 550 nm is chosen despite 
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the fact that its relevant refractive index may vary with 
the ZnO coating method [23,24]. 

We scan the thicknesses of both the photoactive layer 
and the ZnO layer from 5 to 400 nm with a 5 nm step 
for the calculation of the absorbed photon number per 



unit area per unit time, i.e., N p at the photoactive layer. 
This produces data of a 80-by-80 matrix form for the 
three-dimensional graph of N p (T a ,T s ), where T a and T s 
are the thicknesses of the photoactive layer and the op- 
tical spacer, respectively. The absorbed photon number 
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flux per unit time, which is denoted by N p , can be com- 
puted by integration over the standard solar irradiance 
spectrum (AMI. 5 global) [20], as given by 

h 

Np= J j§) [Tt ~ Tb)dX (1) 

where T(A) is the solar irradiance spectrum as a function 
of wavelengths A, X\, and A2 are the lower and upper 
limits of wavelengths for integration, T t and T b are the 
transmissions at the top and bottom of the active layer, 
respectively, and E S (X) is the single photon energy at a 
given wavelength. Then, for a given quantum efficiency 
(QE), the short-circuit current (J sc ) is [25] 

Jsc=N p *QE*e 0 , (2) 

where eo is the elementary charge, eo = 1.6 x 1CT 19 C and 
QE is assumed to be 100%. 



Results and discussion 

Figure 2 shows the calculated J sc as the thicknesses of 
both the photoactive layer and the spacer are scanned. 
The fact that } sc is in proportion to the intensity of the 
field within the photoactive layer reflects that J sc can be 
modulated by thickness modulation of both layers 
through the light intensity redistribution induced by the 
effective optical path-length changes. 

In the absence of a ZnO layer, increasing the photo- 
active layer thickness from 0 to about 80 nm increases 
J sc . The introduction of the spacer layer, however, causes 
this initial increase to occur with a higher slope for lar- 
ger thickness of a ZnO layer (in the ZnO thickness range 
of 0 to 85 nm) as a consequence of the optical field in- 
tensity redistribution at the photoactive layer, due to the 
optical path-length change of the reflected light from 
the electrode. This causes the feature of the improve- 
ment of J sc for the thin photoactive layer by inserting the 
thin ZnO layer, which is similar to [14]. For example, 
decreasing the active layer from 75 to 40 nm produces 





x-position within the device [nm] 

Figure 3 Light intensity profile along the x-direction in the photovoltaic device. Without (a) and with (b) the ZnO optical spacer inserted. 
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Table 1 Complex refractive indices at X = 550 nm for each 
layer in the transfer matrix method calculation 



Number k 



Glass 


1.4516 


0 


ITO 


1.8234 


0.009 


PED0T:PSS 


1.4871 


0.03 


P3HT.PCBM 


2.1944 


0.874 


Optical spacer 


1.4 


0 


Aluminum 


0.96226 


6.68255 



the drop in J sc by only 8.4% upon the insertion of ZnO layer 
of 50-nm thickness. This is due to the ZnO-induced 
change in the interference of incoming photons with the 
photons reflected from the electrode. However, the slope 
at which the J sc increases with increasing the photoactive 
layer thickness becomes small as the ZnO thickness ex- 
ceeds about 100 nm as shown in Figure 2. This indicates 
that for a ZnO layer of thickness ranging from about 
100 to about 200 nm, the increase of photoactive layer 
from 0 to about 100 nm benefits no substantial increase 
in Jsc- 

Figure 3a,b shows the electric field intensity profile 
across the ^-position within the device without and with 
the ZnO layer of 75-nm thickness, respectively. It is vis- 
ible that the insertion of a ZnO optical spacer of a 
proper thickness improves the light intensity at the 
photoactive layer, benefitting absorption of light, there- 
fore, leading to more efficient harvest of light (Table 1). 

The modulation behavior as shown in Figure 2 be- 
comes faded out as the photoactive layer thickness ex- 
ceeds about 300 nm, implying expectation of no further 



significant enhancement of J sc by means of modulation 
of a ZnO spacer thickness. For the study of an optical 
spacer layer of a higher refractive index, we replace the 
ZnO spacer by Ti02 which has the refractive index of 
2.0 at the visible light. We find that the J sc modulates 
more rapidly with respect to the thicknesses in the 
three-dimensional graph (not shown) although the modu- 
lation patterns appear similar to the ZnO case. This sug- 
gests that the Ti0 2 thickness adjustment for the J sc 
improvement (for a given thickness of the photoactive 
layer) requires the finer tuning, leading to more difficulty 
in fabrication of the spacer of such a fine-tuned thickness 
than the case of a ZnO layer. 

It is known that the photoactive layer, i.e., a kind of 
the semiconducting layer, is the most resistant to charge 
transport among the other layers of the device. The 
P3HT:PCBM has the resistivity of about 10 3 Cl ■ cm [26] 
whereas that of a ZnO can be adjusted below 1 Cl ■ cm 
by solution-based method [27] or even below 10~ 3 Cl ■ cm 
by CVD or PVD deposition methods [28]. The electrical 
resistance relevant to the active layer, thereby, must be 
taken into account in the calculation of the device 
efficiency. It is revealed that the device resistance is 
governed by the series resistance consisting of both the 
active layer bulk resistance and the contact resistance 
between the active layer and the Al electrode [15]. It is 
also known that the whole device resistance is approxi- 
mately in proportion to the thickness of the photoactive 
layer, and reducing the resistance leads to higher PCE in 
organic photovoltaic devices. We assume the device is 
annealed in order to decrease the series resistance and 
can thus adopt its relation from [16], which represents 
the linear relationship between the thickness of the 
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active layer {T a ) and the series resistance {R s ) for the 
annealed device, as given by 

R s = 0.06- T a + 3.1. (3) 

This series resistance is used to normalize J sc so that 
we can gain the estimation of the normalized number of 
photons absorbed per unit time per unit cross-sectional 
area of the device, for the device performance, i.e., r\ = 
JsJRs- Figure 4 shows rj (in the unit of mA • cm" 2 • Cl~ 1 • 
cm" 2 ) in a top view of the 3D plot, which is achieved by 
scanning both the active layer thickness from 0 to 
70 nm and the ZnO spacer thickness from 0 to 400 nm. 
The maximum rj that accounts for the most effective 
device performance (r/ = 1.145 x 10 20 mA • cm" 2 • O" 1 • 
cm" , J sc = 10.5 mA • cm" ) occurs at T a = 35 nm and 
T s = 50 nm, while the second maximum occurring at 
T a = 35 nm and T s = 235 nm. 

The characteristics of the monotonic increase in J sc in 
the range of 0 to 70 nm with the spacer thickness 
smaller than 50 nm, shown in Figure 2, are not visible in 
the behavior of the rj, since the rj shows the maximized 
value at the active layer thickness of 35 nm. This indi- 
cates that the series resistance as a function of active 
layer thickness will influence, qualitatively, the depend- 
ence of the device efficiency on thicknesses of both layers 
(the photoactive layer and the spacer layer). We also find 
that minimization of series resistance, e.g., via annealing 
methods, permits not only the relevant voltage to increase 
but also the optimum thickness of the photoactive layer to 
increase, leading to more absorption that can occur. 

Conclusions 

We present, by simulation, the number of photons 
absorbed in the photoactive layer of the BHJ (P3HT: 
PCBM)-based organic photovoltaic device, in the case 
of using an optical spacer inserted between the photo- 
active layer and the electron collecting electrode. 

The adjustment of an optical spacer thickness leads to 
the redistribution of optical intensity of light across the 
photoactive layer, enabling the maximization of the op- 
tical intensity for light absorption maximization, via con- 
trolling optical interference between the forward and 
backward light within the active layer. This interference 
modulation can be achieved by changing the effective 
optical path length of light propagating back and forth 
between the two electrodes, i.e., by changing the thick- 
ness of both the photoactive and spacer layers. We pro- 
vide the 3D plots of the number of photons absorbed 
per unit time per unit area, as a function of thicknesses 
of both layers in the range of 0 to 400 nm. We also in- 
vestigate the effect of using the different index spacer, 
e.g., the higher index layer of TiC>2, and find that more 
delicate tuning of the spacer thickness is required in the 



case of the higher index spacer, for the device efficiency 
optimization. 

In addition, series resistance, i.e., R s , composed of the 
active layer bulk resistance and its contact resistance 
with the electrode is taken into account to provide the 
three-dimensional plot of r; = N p /R s . The results calcu- 
lated show that the series resistance that depends on the 
photoactive layer thickness influences, qualitatively, the 
dependence of PCE on the thicknesses of both layers. 
This reveals that the proper annealing for elimination of 
such thickness-dependent resistance may be vital to in- 
crease not only the relevant voltage but also the optimum 
thickness of the photoactive layer. 

Future works can include the dispersion of refractive 
index and absorption in the constituent layers such as a 
photoactive layer and a spacer for the calculation of the 
device absorption efficiency. 
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